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We investigate the quark backreaction on the Polyakov loop and its impact on the thermodynamics 
of quantum chromodynamics. The dynamics of the gluons generating the Polyakov-loop potential 
is altered by the presence of dynamical quarks. However, this backreaction of the quarks has 
not yet been taken into account in Polyakov-loop extended model studies. In the present work, 
we show within a 2-1-1 flavour Polyakov-quark-meson model that a quark-improved Polyakov-loop 
potential leads to a smoother transition between the low-temperature hadronic phase and the high- 
temperature quark-gluon plasma phase. In particular, we discuss the dependence of our results 
on the remaining uncertainties that are the critical temperature and the parametrisation of the 
Polyakov-loop potential as well as the mass of the a-meson. 
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I. INTRODUCTION 

Understanding the properties of the elementary par- 
ticles in the phase diagram of strongly correlated mat- 
ter is one of the major challenges of present research. 
Its general structure is dictated by the phase transition 
line separating the low-temperature hadronic phase as- 
sociated with quark confinement and chiral symmetry 
breaking and the chirally symmetric high-temperature 
quark-gluon plasma phase. There are several possibilities 
to investigate this transition under different conditions. 
According to the standard model of particle physics and 
cosmology, this transition happened at least once and 
that is in the early Universe at high temperatures and 
small net quark densities (see e.g. Ref. Another 
natural environment where this transition might take 
place, now at high densities and low temperatures, could 
be supernovae and their remnants, compact stars (see 
e.g. Ref. [2]). More detailed information about the phase 
diagram can be obtained by large accelerator facilities 
with high performance detectors, as the Large Hadron 
Collider, the Relativistic Heavy Ion Collider, the Facil- 
ity for Antiproton and Ion Research and the Nuclotron- 
based Ion Collider Facility (see e.g. Ref. [3 ). 

Functional continuum methods are well suited for a 
combined study of the chiral and confining dynamics of 
QCD at finite temperature and density. In recent years, 
much progress in this direction has been made within 
the functional renormalisation group (FRG) approach to 
QCD, see Refs. [HIS], and to Polyakov-loop extended low- 
energy effective models, see Refs. [Sllllj. In fact, the lat- 
ter emerge dynamically from the first-principle QCD-flow 
at low energies: the parameters of these models defined 
below a hadronic mass scale ~ 1 GeV can be extracted 
from a QCD RG-flow starting with the classical QCD ac- 



tion at a given high (perturbative) scale A ^ Aqcd- For 
the Polyakov-loop extended quark-meson model (PQM) 
this has been discussed in Refs. [H El [3 [H] . 

Combining first-principle QCD flows with low-energy 
effective models can be very fruitful. For example, such 
an approach reveals the dynamical connection of con- 
finement and chiral symmetry breaking from an analysis 
of the nonperturbative fixed-point structure of the the- 
ory, see Refs. j9l[10]. More generally, the first-principle 
determination of the input parameters of such models, 
as has been described above, can be used to system- 
atically remove the ambiguities of low-energy effective 
models. This systematic procedure utilises the results 
of the first-principle flows for the chiral properties, see 
Refs. [11[T2HIS], and confinement-deconfinement proper- 
ties, see Refs. [il lTMlQ] . both in Yang-Mills (YM) theory 
and QCD at finite temperature and density. The techni- 
cal details of this systematic embedding of low-energy 
effective models in first-principle QCD is discussed in 
Sec. El 

The embedding discussed above also entails that it is 
even more interesting to study low-energy effective mod- 
els on their own in order to determine the relevant input 
parameters in these models, i.e., those the infrared be- 
haviour of these models is sensitive to. Low-energy ef- 
fective models capture the most important properties of 
strongly interacting matter and are well suited for ba- 
sic investigations and applications, for a recent review 
see Ref. [20] • The most popular ones are the linear sigma 
model, also called quark-meson model in this context [2H - 
132], the Nambu-Jona-Lasinio (NJL) model [33HSZ], and 
the Polyakov-loop extended versions thereof [BHHl HH 1551 - 
[50] . In these models the gauge part is adjusted to lattice 
calculations of Yang-Mills theory [40l|44l|46]. A first step 
towards the full inclusion of the backreactions of quarks 
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on the gauge sector is to estimate the change of the tran- 
sition temperature of the Polyakov-loop potential when 
going from Yang- Mills theory to QCD [6 . This can safely 
be done on the basis of a perturbative estimate of the 
change of Aqcd when going from Yang-Mills theory to 
QCD which also has been used for the chiral dynamics 
[131 HH ■ However, on the mean field level such an anal- 
ysis lacks the full inclusion of glue-matter dynamics to 
both the chiral as well as the confinement-deconfinement 
physics. The related quantum and thermal fluctuations 
typically smoothen the respective phase transitions. In- 
deed, the mean field models show sharper transitions in 
comparison to the full theory, see Ref. |3| . In the present 
work, we will show that the functional renormalisation 
group allows for a comparison of the Yang-Mills poten- 
tial with the glue part of the full effective potential in 
QCD. We extract the relation of both and apply it then 
to improve the Polyakov-loop potential entering effective 
models, such as Polyakov-loop extended Nambu-Jona- 
Lasinio/quark-meson (PNJL/PQM) models. With this 
improved Polyakov-loop potential we shall see that the 
phase transition becomes smoother and is in better agree- 
ment with most recent results from functional QCD stud- 
ies and lattice QCD simulations. 

The paper is organised as follows: in the following sec- 
tion E] we summarise our results from a functional renor- 
malisation group approach which we have already used 
in earlier studies to compute the Polyakov-loop potential. 
In particular, we compare in Sec. |II B] the pure Yang-Mills 
potential with the glue part of the full effective potential 
in QCD. From this comparison, we extract a description 
of how the pure Yang-Mills Polyakov-loop potential en- 
tering PNJL/PQM models has to be modified such that 
the quark backreaction on the gluodynamics generating 
this potential is effectively taken into account. After 
we have given a brief summary of the PQM model in 
Sec. |III[ we employ our results from Sec. |IIB| to amend 
the Polyakov-loop potential entering PNJL/PQM model 
studies. In Sec. |IV[ we then show how this improvement 
affects the temperature dependence of the order parame- 
ters and thermodynamic observables, in particular close 
to the phase transition. To complete our studies we dis- 
cuss the dependence of our results on the remaining free 
parameters. Finally, in Sec. [Vj we give our conclusions 
and outline possible continuations. 



II. FUNCTIONAL FLOWS IN QCD & LOW 
ENERGY EFFECTIVE MODELS 



It has been known since long, see e.g. Refs. [H- 
l51[51 [T5l[?n[nH5i] . that low-energy effective models can 
be systematically related to full QCD within the FRG- 
approach. For Polyakov-loop extended models this fol- 
lows from the Landau gauge approach in Refs. [H [51 155j. 
for the Polyakov gauge see Refs. [T71 [SS]. This set- 
ting has been discussed in details for the PQM model 
in Refs. [51 HI [HI HZ]- While chiral symmetry and its 
dynamical breaking are well described within these mod- 
els, confinement is only included in a statistical manner 
via a phenomenological Polyakov-loop potential lA. This 
potential is fixed to lattice data of the pure Yang-Mills 
system at vanishing chemical potential, hence it is ap- 
proximately the Polyakov loop potential in Yang-Mills 
theory, VVm- More specifically, a Polyakov-loop poten- 
tial is embedded in these models which reproduces the 
temperature dependence of the Polyakov loop and the 
thermodynamics of pure Yang-Mills theory as obtained 
in lattice simulations. As this only fits two observables 
which basically do not change below the phase transi- 
tion, this leaves us with a big parameter space. Results 
for the different potentials have e.g. been discussed in 
Refs. [48l ED] . Also, the coupling of the matter sector to 
the gauge sector is lost in such an approach, see also the 
discussion in Refs. (Hll H [H [SJ . 

In summary this entails that the glue potential of full 
QCD, Fgiuo, encoding the gauge dynamics in the pres- 
ence of matter fields, is replaced by a phenomenological 
Polyakov-loop potential U. The present work is based on 
a qualitative improvement of this approximation towards 
full QCD by embedding Polyakov loop extended models 
within QCD flows. 



A. QCD with functional methods 

In this section we summarise results obtained with 
functional methods for QCD. The present work mostly 
relies on the results obtained in Refs. [H [5j for the QCD 
phase diagram. In Ref. [3] the FRG approach to the 
phase diagram was put forward on the basis of the QCD 
flow equation in the background Landau gauge (Landau- 
DeWitt). It is governed by the flow for the QCD effective 
action rfc[A; </>], 





FIG. 1. Partially hadronised version of the FRG flow for 
QCD. The loops denote the gluon, ghost, quark and hadronic 
contributions, respectively. The crosses mark the RG regula- 
tor term. 



dtTk[A;(l)\ = -TrGa[A;0]ati?A-TrGc[i;0]9ti?c 

-Tr Gq[A- 0] dtRq + ^Tr G^ [A; 0] dtRn , (1) 

depicted in Fig.[l] Here, A is the gluonic background and 
(f> — {a,c,c,q,q, hadrons) are the dynamical fluctua- 
tions with A = A + a. The flow ^ of the effective action 
or free energy only depends on the full field-dependent 
fluctuation propagators G0[^;0], 
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G4A;(l>]{p,q) = 



1 



b,-?), (2) 



and the regulator functions i?^. The regulators R^, de- 
pend on background covariant momenta D = d — i g A 
such that for small momenta it acts as a mass and for 
large momenta it vanishes sufficiently fast; i.e., Rk acts 
as a momentum-dependent additional mass term. The 
present approach has the advantage that hadronic states 
can be included successively within dynamical hadroni- 
sation, see Refs. [51, .52^ i)8, ,59J ■ This simply entails that 
part of the quark and gluonic fluctuations are treated 
separately as they carry hadronic resonances. It is worth 
emphasising that this is done dynamically within the 
QCD-setting so there is no danger of overcounting terms 
present within an effective model setup. The QCD results 
obtained in the present work are computed within the ap- 
proximation to two-flavour QCD detailed in Rcfs. [H [S] 
which explicitly takes into account the lowest mesonic 
multiplet, the cr-meson and the pious, a, tt. 

The FRG-approach to QCD outlined above as well as 
other functional approaches allow us to access the physics 
of the conflnement-deconfinement phase transition in a 
simple way, see Refs. [3J O [1^1 [T51 [TH]. The symmetry 
behind this transition, the center symmetry, is only an 
unbroken symmetry for infinitely heavy quarks. Then, 
QCD reduces to a pure SU(7Vc) gauge theory. The com- 
monly used order parameter is the expectation value of 
the Polyakov-loop operator, 
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Nr. 



Vexpiig J dxoAo{xo,x) 



(3) 



where /3 = 1/T is the inverse temperature, g denotes the 
bare gauge coupling and V stands for path ordering. The 
Polyakov loop is then given by ($ (x)). 

Strictly speaking, the Polyakov loop ($ (x)) is an 
order parameter for center symmetry breaking, see 
e.g. Ref. [60] . However, its logarithm can be associated 
to half of the free energy Fgg of a quark-antiquark pair at 
infinite distance. A center-symmetric confining phase is 
indicated by a vanishing Polyakov loop and implies that 
the free energy of a static quark is infinite. On the other 
hand, the deconfined phase is associated with a finite free 
energy of a static quark and, in turn, a finite Polyakov 
loop. The latter implies that center symmetry is broken 
spontaneously at high temperatures in Yang-Mills the- 
ories. In the presence of (light) dynamical quarks, the 
Polyakov loop is then finite for all temperatures but still 
becomes small in the "confined" low-temperature phase. 
This can be eventually traced back to the fact that quarks 
explicitly break the (global) center symmetry of the the- 
ory 

In addition to the standard Polyakov loop, other order 
parameters for quark confinement have been introduced. 
In this work we shall consider an order parameter which is 
closely related to the Polyakov loop, namely ^[(Ao)]. In 



Polyakov-Landau-DeWitt gauge, it can be indeed shown 
that $[(ylo)] also is an order parameter, see Refs. [T6 l ll7 ) . 
Here, (Aq) is a constant element of the Cartan subalge- 
bra of the gauge group and denotes the ground state of 
the order-parameter potential which is the Polyakov-loop 
potential. Perturbatively, the related effective potential 
l/[(ylo)] has first been computed in Refs. [STl I62j . 

Based on a functional RG approach, first nonpertur- 
bative QCD studies of this potential, including a compu- 
tation of the phase transition temperatures for SU(7Vc) 
Yang-Mills theories, have been presented in Refs. [T5Hl8j . 
In Ref. |19J the computation has been extended to gen- 
eral functional methods, i.e., Dyson-Schwinger equations 
and the 2PI-approach using the thermal propagators ob- 
tained in Refs. [531 [M]. In Refs. [T, the approach has 
been applied to Nf = 2 QCD in the chiral limit. This 
computation within two-flavour QCD includes the full 
backcoupling of the matter sector on the propagators of 
the gauge degrees of freedom via dynamical quark-giuon 
interactions [TH (TSl [H] . In particular, this includes the 
resolution of the full momentum dependence of the ghost 
and gluon propagators [551 - 155] . Recently, the nonper- 
turbative Polyakov-loop potential has also been studied 
using various other approaches [7DH75]. 

In the present work, we refrain from repeating the de- 
tails of the above works but only refer to it whenever it is 
required. For general QCD-related introductions and re- 
views to our functional RG approach, we refer the reader 
to Refs. |5l|23l52H5l[7Il[I8l. 



B. 



Prom QCD to Polyakov-loop extended chiral 
models 



For temperatures below the chiral and conflnement- 
deconfinement phase transitions the glue dynamics de- 
couples from the matter dynamics. In the Landau gauge 
this physical decoupling is realised simply by a mass gap 
in the gluon propagator, that is p^Ga{p^) — > for small 
momenta, see e.g. Ref. 'W and references therein. As the 
ghost-matter coupling is mediated by the gluon this de- 
coupling extends to the full ghost-gluon dynamics. Hence 
the pure glue sector of the theory decouples as expected, 
and the QCD flow equation reduces to 

dtTklA; 0] = -Tr GJA; OtRg + ^Tr Gh[A; </)] OtRn , 

(4) 

which is depicted in Fig. [2] The background gluon fleld 
as well as the ghost and gluon fluctuation flelds A and 
c, c are simply spectators in Eq. Q. If setting the specta- 
tor fields to zero the flow equation Q reduces to that of 
generic low-energy models of QCD. The respective mod- 
els are then singled out by specifying the hadronic con- 
tent of the flow. In Eq. Q, the presence of a quark loop 
should not be confused with the presence of quarks as ob- 
servable (asymptotic) states. The presence of this loop 
only states that quark loops are still present "internally" 
and required to give a microscopic description of the in- 
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dtTM = 



FIG. 2. FRG flow for the matter sector of QCD. The loops de- 
note the quark and hadronic contributions, respectively. The 
crosses mark the RG regulator term. 



teraction of the hadronic degrees of freedom. In the same 
spirit, gliion loops also play a role in the low-energy sector 
even though they do not represent observable (asymp- 
totic) states in this regime as well. However, gluon loops 
govern the dynamics of the confinement order-parameter 
potential which in turn determines the background gluon 
field entering Eq. Q, see also our discussion below. 

The most prominent example are NJL-type models. In 
these models the hadronic loops are missing completely 
and the full matter dynamics is described by the quark 
loop. This is in correspondence to a description of the 
low-energy matter dynamics of QCD purely in terms of 
quark correlation functions. This only captures the cor- 
rect QCD-dynamics if the hadronic spectrum is taken 
into account via resonances in the scattering amplitudes 
of quarks, e.g. in the quark-quark scattering kernel. The 
latter is directly related to the four point function T*^^^ of 
two quark-antiquark pairs. It is worth noting in this con- 
text that the flow equation Q includes vertex and prop- 
agator corrections to the theory and hence naturally in- 
corporates the so-called nonlocal NJL models [56 1 179 1 180]. 
For a detailed discussion of the presence and emergence 
of nonlocal fermionic vertices in RG flows, we refer the 
reader to Ref. [51) . 

Even though this incorporates in principle the full dy- 
namics, hadronic contributions are more easily taken into 
account via (dynamical) hadronisation, see Refs. [151 [HI 
[52l l58l [59] . This technique describes the following physi- 
cal mechanism: when lowering the RG-scale in the QCD 
flow equation ([T]) , Fig. [ij strongly-bound quark correla- 
tions e.g. quark-antiquark (mesons), diquark, and three- 
quark (baryons) correlations, will be dynamically created 
by the flow. These correlations can be parameterised in 
terms of mesonic, diquark and hadronic operators respec- 
tively. Note that the diquark operators do not describe 
asymptotic states but intermediate correlations. The re- 
spective coupling of the hadronic operators is dynami- 
cally enhanced by the quark-gluon fluctuations at higher 
scales, while it takes over the dynamics from the quark- 
gluon sector at lower scales, hence the name dynamical 
hadronisation. This fluctuation-triggered hadronisation 
within the QCD-flow clearly avoids any double-counting 
problem: the hadrons are created dynamically, the re- 
spective couplings are created from the QCD input which 
only depends on a single input parameter, the strong cou- 
pling as- 



Quantum and thermal fluctuations are self-consistently 
taken into account. Still, at low scales, say for fc < 1 GeV, 
the pure glue sector decouples and we are left with 
Eq. Q. Note that within such a procedure the initial 
effective action Fa at A « 1 GeV serves as the classical 
action of the quark-hadronic low-energy model and its 
coupling constants are directly derived from QCD. 

In the present work we consider such a quark-hadronic 
model at finite temperature and vanishing density: close 
to the phase boundary between the quark-gluon plasma 
phase and the hadronic phase and at not too large chem- 
ical potential mesonic degrees of freedom, in particular 
the pion and d-fluctuations, become important. Here, 
we assume that the associated bosonic fields can be 
described as composites of fermions and do not carry 
an internal charge, such as colour or flavour: a ^ qq 
and if ^ 975 tq. A priori, these fields are just auxiliary 
fields introduced by means of a Hubbard-Stratonovich 
transformation of a purely fermionic theory. Since these 
composite objects have considerable overlap with the full 
meson operators, they are well suited to describe the me- 
son dynamics in a first approximation. Keeping this in 
mind, the (effective) action of QCD can be expanded 
in powers of these mesonic operators, at least for small 
momentum scales. Momentum dependencies and correc- 
tions of higher order can conveniently and systematically 
be taken into account, e.g., by means of a derivative ex- 
pansion of the effective action, see Ref. [S3] for a review. 
In the fully coupled QCD RG-flow, the strength and the 
momentum dependence of the various couplings can then 
be derived from the microscopic QCD action. In fact, the 
loops depicted in Fig. [l] are intrinsically coupled to each 
other. In particular, there are contributions from the 
matter sector to the diagrams for the gluon propagator, 
i.e., the quark part of the gluonic vacuum polarisation, 
see Fig. [4] 

Hence we only take into account mesonic correlations 
in terms of effective operators. Baryonic fluctuations are 
considered as subleading at vanishing density and are 
partially taken into account in the quark correlations. 
Moreover, the gluonic background is determined by the 
constant solution of the QCD equations of motion. Such 
a solution can always be rotated in the Cartan, 



(5) 



with constant {Aq'^). The minima and maxima of the 
potential are accessed for (Aq) — 0, see Refs. [^ irMTO] . 
Then the Polyakov loop reads 



mo] 



2 cos ( -pgAo 



with Ao = (Al) (6) 



In summary this leads to quark-meson flows Q in con- 
stant temporal gauge field backgrounds which define the 
PQM model [5] , including its matter quantum and ther- 
mal fluctuations. We note that for the NJL-type model 
this deflnes the PNJL model including its quantum and 
thermal fluctuations. The generalised version is the 
Polyakov-loop extended quark-hadron (PQH) model. 
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FIG. 3. Functional flow for the effective action. Lines with 
fiUed circles denote fully dressed field-dependent propagators 
([2|. Crossed circles denote the regulator insertion dtRk- 



Even though the glue-dynamics decouples for low tem- 
peratures, the gluonic background A — (A) does play 
an important role. In the confining phase it screens the 
propagation of the quarks which confines them even sta- 
tistically, leave aside the confining potential. The ex- 
pectation value (A) is given by the solution of the QCD 
equations of motion which also requires the first part of 
the QCD flow, that is 



dtrk[A;<j>] = ^TrGaiA; 



\dtRA-TrGc[A;4,]dtRc. (7) 



depicted in Fig. |3] The fluctuation is evaluated at the 
equation of motion (EoM) , (f) = (f>: the gluonic fluctuation 
background is vanishing, a = 0, whereas the mesonic 
background is given by cr = a and tt = leading to the 
running quark masses in the quark vacuum polarisation 
depicted in Fig. |4j This leads us finally to 



dtVsiue[Ao] = —dtTk[Ao;4>] 



(8) 



for constant backgrounds Aq and spatial volume V. 
Structurally, Eq. Q resembles the pure Yang-Mills flow 
equation, 

dtVyulAo] = -^dtTYMA^o; a = 0,C = 0,C = 0], (9) 

which also has the pictorial form of Fig. [3] The ghost 
and gluon propagators in Eq. ([8|, however, are those of 
QCD. In particular, the flow of the gluon propagator re- 
ceives contributions from matter loops, e.g. the quark 
contribution to the vacuum polarisation, see Fig. |4j 

It has been already mentioned at the beginning of this 
section that in the standard approach to Polyakov-loop 
extended models the glue potential Vgiuc in Eq. (jsl) is ap- 
proximated by an approximation to its Yang-MiUs ana- 
logue, Vym in Eq. To bring these studies closer to 
(full) QCD, however, this Polyakov-loop potential effec- 
tively generated by pure gluodynamics has to be replaced 
by the QCD glue potential, i.e., by the contribution stem- 
ming from the gauge degrees of freedom in the presence 
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FIG. 4. Quark polarisation contribution to the gluon propa- 
gator representing a contribution to the matter backcoupling. 
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Yang-Mills 
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PgAo / 2k 

FIG. 5. Comparison of the SU(3) Yang-Mills and glue ef- 
fective potentials as functions of the background gauge field 
Aq = (Aq) for various reduced temperatures, see the text 
for the definition of the reduced temperatures. The form of 
the potentials is very similar, however, the temperature scale 
changes. Only the glue part of the full effective potential is 
shown here and compared to the Yang-Mills potential. 



of dynamical quarks. It is therefore beneficial to amend 
these model calculations by utilising the available infor- 
mation on this glue part of the potential. 

In the present work, we only consider Polyakov-loop 
potentials which have been computed by considering 
two quark flavours with vanishing current masses. The 
strange quark has been ignored. We expect that the 
corrections resulting from the inclusion of a strange 
quark are subleading in our present studies. Following 
Refs. [H [TBHTO] . we have computed the nonperturba- 
tive Polyakov-loop potential using the above mentioned 
functional RG approach taking into account the back- 
reaction of the quark degrees of freedom on the gluon 
propagators [Ij. 

In Fig. [5] we compare the effective potential of SU(3) 
Yang-Mills theory as obtained in Refs. 16, 18, 19] with 
the Polyakov-loop potential of the glue sector of our 
present study. The latter includes also the quark part 
of the gluonic vacuum polarisation but does not include 
the fermionic part of the full potential. The related order- 
parameters $[(v4o)] derive by inserting the temperature- 
dependent gauge field value at the minimum into Eq. (|6| . 
The order parameter shows a first-order phase transition 
for both, the glue potential and the Yang-Mills poten- 
tial. It is displayed for Yang-Mills theory in Fig. |6] We 
observe immediately in Fig. [s] that the temperature scale 
of the various Yang-Mills potentials and the glue poten- 
tials differ. The differences are induced by the matter 
fluctuations altering the propagators of the gauge fields. 
This has already been anticipated in Ref. [6], where a 
phenomenological hard-thermal-loop (HTL) estimate led 
to T^.^^ = 208 MeV for two quark flavours. 
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FIG. 6. Yang-Mills order-parameter <E>[{j4o)] as a function of 
the normalised temperature T/Tc with Tc = 276MeV. 



Although the temperature scales are different, the 
shape of the order-parameter potential V is not, see 
Fig. [5] We can exploit this observation to estimate how 
the temperature of a given pure Yang-Mills Polyakov- 
loop potential has to be modified to be closer to the QCD 
potential. To this end, we need to define a measure for 
the comparison of the potentials. We use 



Moreover, we introduce the reduced temperatures 



glue 



^ cr 



(10) 



(11) 



where rs.'"'= = 203 MeV. Note that this is not the criti- 
cal temperature obtained from the full effective potential 
including pure quark loops but only from the contribu- 
tion stemming from the glue part of the potential for two 
massless quark fiavours [J]. We add that the absolute 
scale in Ref. ^ was not computed in a chiral extrapo- 
lation of the theory with physical quark masses. Hence, 
in the applications to the PQM model we will leave the 
glue critical temperature as a free parameter, maximally 
ranging from 



180 MeV < rs'"° < 270 MeV . 



(12) 



The upper limit in Eq. ( |12[ ) is the critical temperature 
of Yang-Mills theory, the lower limit is the estimate in 
Ref. H], see also Table [l] For T™, on the other hand. 



TABLE I. Critical temperature of the glue effective potential 
Tc^r"" for ^ ! massless flavours and in case of an additional 
massive quark flavour with the current strange quark mass of 
95 MeV [81], according to Ref. [6]. 
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FIG. 7. Relation between the two temperature scales of pure 
Yang-Mills theory and the glue part of our present study in- 
cluding the backreaction of the quarks on the gauge flelds. 
The solid (green) line is our fit to the numerical data (black 
crosses) in the range around the transition temperatures. A 
clear deviation from linear behaviour with slope one (dashed 
red line), i.e., data from pure Yang-Mills theory, can be recog- 
nised. 



we find T™ = 276 MeV [H [H HH] which is in quan- 
titative agreement with the results from lattice studies. 
The comparison of the potentials then yields the trans- 
lation of the two temperature scales: given a reduced 
QCD-temperature igiuo, the related Yang-Mills temper- 
ature ^YM is that which minimises Eq. (10 1. This leads 
to 



iYM(^gluc) ~ 0.57tglue : 



(13) 



which, together with the absolute temperature scale in 
Eq. ( 12 ) serves as an important input for model studies. 
This relation is displayed in Fig. [t] Equation ( 13 ) also 
provides the map from the temperature-dependence of 
the Yang-Mills order-parameter to that derived from the 
glue potential by using tYM(igiuc) in Fig. [6] The respec- 
tive order-parameter fiattens, but still shows a first-order 
transition. This procedure also extends to finite chemical 
potential e.g. with the chemical potential dependence in 
Ref. |6]. A refinement of this estimate based on full QCD 
flows, which obeys the Silver-Blaze property for vanishing 
temperature, is derived in Ref. [11]. Note also that the 
approximation ( 13 1 holds only for small and moderate 



temperatures, as at high scales the slope of Eq. ( 13 1 sat- 
urates, where one reaches the perturbative limit and the 
potentials reach their asymptotic form. In fact, we would 
naively expect that the results for the reduced tempera- 
tures agree at (very) high temperatures where the quark 
degrees of freedom are parametrically suppressed. 



In the following we show how our prescription ( 13 ) 
can be used to amend PNJL/PQM model studies. More 
precisely, we use this prescription to improve the poten- 
tials for the Polyakov loop entering these models. Al- 
ternatively to our procedure outlined in the subsequent 



7 



sections, we could also directly use the Polyakov-loop po- 
tentials obtained from our functional RG approach in a 
PNJL/PQM model study to get even closer to the full 
QCD study in Refs. [H [51 [SS]. For first studies in this 
direction, we refer the reader to Refs. ^ JJ3 where the 
ground state of the nonperturbative Polyakov-loop po- 
tential has been used to study the relation of confinement 
and chiral symmetry breaking. 



TABLE II. Values of constants to which the parameters of the 
mesonic potential are adjusted to, according to Ref. ^1] and 
value of the constituent quark mass of the light (up and down) 
quarks that we use to fix the quark-meson Yukawa coupling 
in Eq. ([l8|. 



Constant 


U 


/k 


rrin rriK m^/ ttio- 


mi 


Value [MeV] 


92 


110 


138 495 548 958 400-600 


300 



III. POLYAKOV-QUARK-MESON MODEL 

A. PQM model with 2+1 quark flavours 

The particle content of the PQM model with 2-1-1 
flavours are the constituent quarks minimally coupled 
to gauge fields, and coupled to mesons via a Yukawa- 
type term. The coupling to the mesons generates the 
masses of the quarks because the mesonic potential con- 
tains spontaneous and explicit chiral symmetry break- 
ing. The action of the PQM model with 2-1-1 flavours 
includes kinetic terms for the dynamical low-energy de- 
grees of freedom, quarks and mesons. The quarks are 
minimally coupled to a background gauge field (Aq) as 
discussed in Sec. HIl which can be rewritten in terms of a 
coupling to the Polyakov loop background, <i>[(^o)]- The 
kinetic part of the model is complemented by the effec- 
tive potential as a function of the order parameters and 
thermodynamic control parameters. In the present work 
we restrict ourselves to isospin symmetric matter, and do 
not distinguish between the up and down quark sectors. 
Hence, the order parameters of chiral symmetry are the 
light or nonstrange condensate a and the strange chi- 
ral condensate CTs and the pion expectation values vanish 
7? = 0, e.g. Refs. (13 ED- Additionally the pure glue 
sector in Eq. ([I]) , Fig. l] leads to an effective potential 
of the Polyakov loop $"(]4o)]. In summary, the effective 
potential of the model reads 

n (a, a„ $; T, /i/) = [/ (a, a,) + U T) + 

+ r!qq(a,a„$,$; T,fif ) , (14) 

where $ is the conjugate of The mesonic order- 
parameter potential for chiral symmetry breaking can be 
written as follows (13 [32 



U [a, CTs) 



1 



1 



+ - (2Ai + A2) CT* + - (2Ai + 2A2) (7, 



2^2 



— ha — h,a^ 



(15) 



Hence, the mesonic sector has six parameters that are 
the coupling constants m^, Ai, A2 and c and the ex- 
plicit chiral symmetry breaking terms h and h^. They 
are adjusted to the pion and kaon decay constants /^r 
and /k and the meson masses of the scalar and pseu- 
doscalar octet, namely TOtt, "t-k, to^ + 'm'ff 



The mass of the sigma meson is not exactly known. 
Within our model the sigma meson is identified with the 
experimentally measured resonance /o(500) with mass 
TOf„ = (400- 550) MeV [HI]. We wiU discuss the de- 
pendence of our results on this uncertainty. Once this 
set of masses and decay constants is given, the model 
parameters are defined. For their explicit derivation we 
refer to Refs. (13 [21]. The values of the constants we use 
to calculate these parameters are listed in Table [llj 
The last term of Eq. ( 14 ) represents the constituent 



quark sector and derives from the QCD flow equation ([I]), 
Fig. [1] within a 1-loop approximation to the quark loop. 
It includes the coupling to the Polyakov-loop variable 
$[(^0)] and the mesons 



f]qq(a,as,$,$;T,^/) = -2T ^ 



(2^ 



/=u,d,s ' 



X |ln 1 + 3 ( 



In 



1 + 3 + ^e-t-^^+^^^/^j X 
The dispersion relation of the quarks, 



Ef = Jk^ 



(16) 



(17) 



couples the chiral condensates to the Polyakov-loop vari- 
ables since the constituent light and strange quark masses 
are, respectively. 



mi 



and 



g 

T2 



(18) 



The Yukawa coupling, g, in Eq. ( |18[ ) between quarks 
and (pseudo-)scalar mesons is fixed by choosing the con- 
stituent mass of the light quarks to be m\ — 300 MeV, 
which results in mg ~ 417 MeV for the constituent 
strange quarks. 

One can combine nonstrange or light and strange con- 
densate to what is called the subtracted condensate 



M.s 



(19) 



T=0 
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This quantity is better accessible in lattice simulations 
than the condensates themselves and hence will be used 
for the comparison of the behaviour of the chiral sector 
in the PQM model with lattice calculations. 

It can be easily inferred from the FRG-setting dis- 
cussed in Sec. |ll] that the effective potential detailed 
above has to change under a variation of the physical 
UV-scale A of the low-energy model. This follows from 
the fact that A also plays the role of an IR scale, up to 
which QCD fluctuations are integrated out, see Sec.|llj In 
an FRG-setting the A-dependence can be directly com- 
puted from the flow and follows from the invariance of 
the effective action at vanishing cutoff scale under a vari- 
ation of A, see the reviews [3 EH EHH [72l [TOl [82H84]. 
In the model setup this well-known fact is related to the 
so-called vacuum (or sea) contribution to the effective 
potential, see Ref. [85]. Its origin in the FRG-setting 
makes clear that it already includes part of the quantum 
fluctuations related to the UV-scale A. Consequently we 
expect significant contributions in particular for larger 
chemical potential, large mass scales and higher temper- 
ature. The impact on the effective potential in the PQM 
model is discussed in more detail in Ref. [IT]. In the 
present work it is a subleading effect as we only consider 
vanishing densities and adopt low values for the mass of 
the cr-meson. 



B. Polyakov-loop potential 

The Polyakov-loop potential U ($, T) simply is the 
pure glue potential V^iue[(^o)] defined in Eq. ([s]). How- 
ever, in the standard approach to Polyakovdoop ex- 
tended models this potential is usually modelled in 
the following way: The functional form for the poten- 
tial contains all terms invariant under transforma- 
tions up to quartic order in the Polyakov loop (see 
e.g. Refs. [SnmOlillllllS]). One possible parametri- 
sation of the Polyakov-loop potential is the polynomial 
parametrisation of Refs. [40l|44]. Reference [40] used a 
different definition of the coefficients than Ref. [33]. A 
simple calculation allows the translation from one set of 
coefficients to the other, see Ref. l50l. 



2^4 



2 6 ^ 



+ T (20) 



with the temperature-dependent coefficient &2 defined as 

as 



i + < {1 + ty {i + ty 



(21) 



TABLE III. Parameters of the gauge potential parametrisa- 
tions for fits to the lattice Yang-Mills simulations [86l |87] . 





ao 


m 


02 


as 


63 


&4 


Poly-I 40 


1.53 


0.96 


-2.3 


-2.85 


13.34 


14.88 


Poly-lf 41, 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 




Ao 


Ai 


A2 


Ba 






Log [ig 


3.51 


-2.47 


15.2 


-1.75 







A{t) 



rp4 

X In 



-t- B (t) X 



1 - 6$$ + 4 ($^ 



3($$)^ 



where both coefficients are temperature dependent, 

Ai A2 



A (<) = Ao + 



1+i {l + tf 



B{t) = 



{1 + ty 



(22) 

(23) 
(24) 



The parameters of the model potential are then deter- 
mined with the help of Yang-Mills lattice data [SHUHZ] for 
the Polyakov-loop expectation value ($), that is the loca- 
tion of the absolute minimum, and the thermodynamics, 
that is the value of the potential at the minimum. They 
are listed in Table Hill 

To compare the different parametrisations and param- 
eter sets and for the later discussion of their impact we 
show in Figs.[8|and|9]the location of the minima of the po- 
tential and the normalised trace anomaly. We see that, 
first, the logarithmic parametrisation is in best agree- 
ment with the lattice data, second, the Poly-I potential 
features the strongest phase transition, third, the Poly-II 
potential leads to a relatively weak first-order transition 
and shows an offset of the peak of the trace anomaly com- 
pared to the other two potentials and the lattice data. 
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Another possible parametrisation of the Polyakov-loop 
potential of Ref. [46] contains a logarithmic term that re- 
stricts the Polyakov-loop expectation value to be smaller 
than one, 



FIG. 8. Polyakov loop as a function of the reduced temper- 
ature in the pure-gauge sector for the different parametri- 
sations and parameter sets of Table |III[ compared to corre- 
sponding lattice results taken from Ref. |§7) . 
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reduced temperature tyy 

FIG. 9. Scaled trace anomaly as a function of the re- 
duced temperature in the pure-gauge sector for the different 
parametrisations and parameter sets of Table |III[ compared 
to corresponding lattice results taken from Ref. 



The different description of the strength of the phase 
transition shows a significant impact on the surface ten- 
sion in the PQM model at finite density [50l I89j. 

The construction of the Polyakov-loop potential 
U entails that it models the pure gauge potential 
Uym (^ym) /Tym in terms of the variable ($). It de- 
scribes successfully the first order phase transition of 
SU(3) Yang-Mills theory. The critical temperature of 
the Polyakov-loop potential entering the reduced temper- 
ature is accordingly that of pure gauge theory, T^,^ = 
270 MeV. 

We conclude that the use of the model potentials 
U incorporates several approximations and deficiencies: 
Firstly, a Yang-Mills potential is used instead of the glue 
potential. Secondly, it only models the Polyakov-loop po- 
tential in Yang-Mills theory by using the location of its 
absolute minimum and the value at the minimum. Below 
Tc this basically provides no information about the po- 
tential. Thirdly, the potential has to be known as a func- 
tion of <I>[(Ao)] as the coupling to the matter sector is de- 
scribed in this variable. Strictly speaking, ^[(ylo)] « (<&) 
is only valid in the Gaufiian approximation, in general 
we have $[(Ao)] > ($) up to renormalisation issues, see 
Refs. [Ullij. 

In the present work we resolve the first approxima- 
tion by using the relation between pure gauge and glue 
effective potential ([Is]). This allows us to convert the 
Yang-Mills potential Uym for the Polyakov loop to a glue 
potential of full QCD, 



glue 
2^4 



luc) — 



Uym 



;*,$,iYM(igiuc)) , (25) 



with tYM(^giuc) in Eq. (13 1. The absolute temperature 



scale is then determined by the value of the transition 
temperature ( [T2| . The resolution of the other two ap- 
proximations will be discussed elsewhere. 



IV. RESULTS AND DISCUSSION 

Given the temperature T and the quark chemical po- 
tentials fi f , the effective potential ( 14 ) is given as a func- 
tion of the four order-parameters a, as, $ and 4>. In 
equilibrium, the expectation values of the order parame- 
ters are given by extremising the effective potential 



da. 







(26) 



For vanishing density the EoMs in Eq. (26), with the 
constraint <!> = $, are given by a minimum of the effective 
potential. In contradistinction, for nonvanishing density, 
with <I> 7^ $ the solution of Eqs. (26) is only a saddle 



point, see the discussions in Refs. f6t i50i i90j . 

In this work we compare the solutions of the equations 



of motion ( 26 ) using for the Polyakov-loop potential the 
pure YM potential and the enhancement to a glue po- 



tential (25) by the relation (13) 



The lattice results we compare our results to are the 
continuum extrapolations of Refs. (HH El] and those of 
the HotQCD collaboration using the HISQ action and 
temporal lattice extent N^. = 12 (dots) and N^. = 8 
(squares) of Refs. [131 Hi], both with physical quark 
masses. 

To investigate the improvement of the Polyakov-loop 
potential by Eqs. ( 25 ) and ( 13 ) we first show results with 



the combination of glue critical temperature, mass of the 
sigma meson and parametrisation of the Polyakov-loop 
potential that reproduces best the lattice results. After- 
wards, we discuss the dependence of our results on the 
uncertainties of these parameters. 



A. Improved Polyakov-loop potential 

To compare the results of the PQM model using 
the Yang-Mills Polyakov-loop potential and the quark- 
improved Polyakov-loop potential we first have to choose 
a canonical setting for the free parameters and uncertain- 
ties that are the critical temperature and the parametri- 
sation of the Polyakov-loop potential as well as the mass 
of the (T-meson. 

The logarithmic parametrisation of the Polyakov-loop 
potential ( 22 ) with the parameters of Ref. |46| is in best 
agreement with the expectation value of the Polyakov 
loop and the equation of state of Yang-Mills theory as 
can be seen in Figs. |8] and [9j Therefore, we use this 
parametrisation for the following discussion. Further- 
more, we take a medium value for the mass of the sigma 
meson of mo- = 500 MeV. To adjust the critical tempera- 
ture of the glue effective potential in full QCD we want to 
bring our model in agreement with lattice results of the 
normalised pressure p/T* since the pressure is directly 
related to basic quantities in lattice simulations (the free 
energy) and mo del studies (the minimum of the poten- 
tial). Figureliolshows that with Tg'"= = 210 MeV we find 
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FIG. 10. Scaled pressure "pjT^ as a function of temperature at 
jif — 0. The red dashed hne is the resuH using as Polyakov- 
loop potential the pure Yang-Mills potential, the full green 
line is the case for the improvement to the glue potential. 
The orange band is the fit to the continuum extrapolation of 
the lattice results of Ref. [92] . 



a normalised pressure p/T^ that is in very good agree- 
ment with the continuum extrapolated lattice results of 
Ref. [92]. 

Now that wc have fixed all uncertainties we can inves- 
tigate how the results evolve when we amend the pure 
Yang-Mills Polyakov-loop potential to a glue potential 
(25) by the relation (fTsl). In Fig. 10 we see that the en- 



hancement of the Polyakov-loop potential smoothes the 
phase transition significantly bringing the effective model 
in close agreement with lattice and functional calcula- 
tions. This agreement holds in the low as well as in the 
high temperature phase. To achieve this agreement in 
the phase were chiral symmetry is broken we took into 
account the contribution of thermal pions. 

The critical temperature of the glue potential Tf^^° = 
210 MeV for which we find the best agreement with 
the lattice calculation is larger than the estimate for 
2-1-1 quark flavours of Ref. |BJ. But this estimate of 
rpgue (2 + 1) = 182 MeV is only a lower limit as outhned 
in Sec. [n] when discussing Eq. ( |l2| ). We find compara- 
ble curves for the combinations (mo- = 400 MeV, T£^'^ = 
180 MeV) and {m„ = 600 MeV, Tg""" = 250 MeV). We 
will discuss the impact of the critical temperature of the 
Polyakov-loop potential and of the mass of the cr-meson 
in detail in the following sections. 

To complete the discussion of the pressure shown in 
Fig. [To] let us note that the result of the PQM model with 
the quark-improved Polyakov-loop potential still shows a 
transition that is marginally steeper than that on the 
lattice. A quantitative probe of this slight deviation is 
the trace anomaly or interaction measure which tests not 
only the free energy or minimum of the potential but con- 
tains as well information about their temperature deriva- 
tive. In Fig. [11] we compare results for the normalised 
trace anomaly (e — 3p) /T^. Even though the agreement 



FIG. 11. Scaled trace anomaly (e — 3p) /T'* as function of 
temperature at = 0. The red dashed line is the result us- 
ing as Polyakov-loop potential the pure Yang-Mills potential, 
the full green line when the improvement to the glue potential 
by the adjustment from functional calculations is taken into 
account. The black band is the fit to the continuum extrap- 
olation of the lattice results of Ref. [92] and the data points 
lattice results of the HotQCD collaboration [94] . 



of the normalised pressure of the effective model using 
the improved glue potential and the lattice calculation is 
remarkably close the slight differences get more visible in 
the trace anomaly. The amplitude is in good agreement 
with lattice results, nevertheless the transition in the ef- 
fective model is still steeper than in lattice calculations. 
However, the calculations done here within a basic effec- 
tive model neglect contributions from correlations. Tak- 
ing those into account smoothes the phase transition, as 
is shown e.g. in Refs. [5] US] which will lead to a better 
agreement. 

Anyhow, one sees in the interaction measure that the 
improved glue potential smoothes the transition com- 
pared to the standard calculation with the Yang-Mills 
Polyakov-loop potential. So the emergence of the quark 
degrees of freedom happens in a larger temperature inter- 
val and its amplitude is in better agreement with lattice 
results. 

After this analysis of thermodynamic quantities, let 
us now analyse the evolution of the order parameters 
that are the Polyakov loop for deconfinement and the 
subtracted chiral condensate ( 19 1 for chiral symmetry 
restoration. 



The biggest discrepancy between Polyakov-loop ex- 
tended effective models and lattice results for full QCD 
is seen in the Polyakov-loop expectation value [Hj. The 
lattice data shows a smoother transition with significant 
smaller values. This can be seen in Fig. [12] and holds as 
well when the quark-improved Polyakov-loop potential is 
applied. Including contributions of fluctuations to the 
model will certainly reduce this large offset to the lattice 
data. However, a part of the discrepancy originates in 
the inherent approximations which are still present: the 
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FIG. 12. The Polyakov-loop expectation value $ as function 
of temperature at ^/ = 0. The red dashed line is the result 
using as the Polyakov-loop potential the pure Yang-Mills po- 
tential, the full green line when the improvement to the glue 
potential by the adjustment from functional calculations is 
taken into account. The grey band is the continuum extrap- 
olation of the lattice results of Ref. [91] and the data points 
lattice results of Ref. . 



FIG. 13. The subtracted chiral condensate Ai,s as a func- 
tion of temperature at — 0. The red dashed line is the 
result using as Polyakov-loop potential the pure Yang-Mills 
potential, the full green line when the improvement to the 
glue potential by the adjustment from functional calculations 
is adopted. The grey band is the continuum extrapolation 
of the lattice results of Ref. [5T] and the data points lattice 
results of Ref. [93| , see also Ref. [95] . 



derivation of the PQM model from QCD entails that the 
Polyakov-loop variable in the quark loop is $ [(^o)] and 
not ($) as used in the Polyakov-loop-model potentials 14. 
This mismatch can be resolved by using the QCD glue 
potential Vgiuc as defined in Eq. (j8|. Then the model is 
fully consistent with its QCD counterpart defined by the 
flow equation ([T]) . The Polyakov loop observable accessi- 
ble in such a continuum approach to full QCD is $ [(^o)] 
and not ($) computed on the lattice. The former observ- 
able has been computed in Ref. [1] within QCD-flows for 
two-flavour QCD in the chiral limit. It shows a rather 
smooth transition but cannot be directly compared with 
present lattice data. On the lattice, the chiral limit is not 
yet accessible. In any case, the two observables are differ- 
ent and mapping them into each other is an interesting 
open issue. 



TABLE IV. Pseudocritical temperatures for the crossover 
phase transition at =0. They are determined by the peaks 
in the temperature derivatives of the subtracted condensate 
Ai,s. 



One sees in Fig 12 that by applying the improved glue 
potential, the Polyakov loop is shifted differently in the 
two phases, as anticipated in Figs. [5] and [t] and Eq. (13 1. 
So, the Polyakov-loop variables are shifted to higher ex- 
pectation values in the confined phase and to lower val- 
ues in the deconfined phase. This shift in the respective 
phases is what leads to a overall smoother transition, not 
only of the Polyakov-loop expectation value. The evolu- 
tion of the Polyakov loop when the quark-improved po- 
tential is considered has the advantage that the Polyakov- 
loop variables are shifted to a value of ~ 0.5 at the tran- 
sition which is a reasonable criterion to define the tran- 
sition between the confined phase and the quark gluon 
plasma in case of a crossover. 

The adjustment of the Polyakov-loop potential from 
the pure gauge potential to the glue potential in full QCD 
does not only affect the Polyakov-loop expectation value 





Wym 
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Budapest [91] 


lattice HotQCD 
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Tc [MeV] 


168 


158 
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159 & 163 



but also quarks and mesonic degrees of freedom because 
of the nontrivial coupling to the gauge field. Therefore, 
we show in Fig. [TB] how the subtracted chiral condensate 
evolves with increasing temperature for vanishing den- 
sity. The result of the PQM model using as Polyakov- 
loop potential the pure gauge potential and the improve- 
ment to the glue potential from functional calculations 
show both a steeper decrease in the transition region 
than the lattice calculations. Nevertheless, the adjust- 
ment of the gauge potential to the glue potential leads 
to a smoother chiral transition and to an improvement 
towards the result of lattice calculations. 

In the previous and following figures the abscissae 
are in units of the reduced temperature of full QCD 
t = (T—Tc) /Tc- To have a unique criterion for lattice and 
model calculations we determine the pseudocritical tem- 
perature Tc of the crossover transition by the peak of the 
susceptibility of the subtracted chiral condensate. Table 
jl Vj summarises the hereby defined pseudocritical temper- 
atures for the lattice calculations and the model calcula- 
tions with the Yang-Mills and improved glue Polyakov- 
loop potential. We see that applying the quark-improved 
Polyakov-loop potential leads to a reduction of the pseu- 
docritical temperature. With Tc = 158 MeV we get a 
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FIG. 14. Scaled pressure vjT^ as a function of temperature 
at /i/ = with the adjustment to the glue potential for four 
different glue critical temperatures. The brown band is the 
fit to the continuum extrapolation of the lattice results of 
Ref. [H]. 



FIG. 15. Scaled trace anomaly (e — 3p) /T* as function of 
temperature at fj,f =0 with the adjustment to the glue po- 
tential for four different glue critical temperatures. The black 
band is the fit to the continuum extrapolation of the lattice 
results of Ref. [92] and the data points lattice results of the 
HotQCD collaboration [94|. 



result that is in very good agreement with the best lat- 
tice calculations. That is another important ingredient 
to find the nice agreement of our results for the equation 
of state and order parameters. 

We conclude that adjusting the gauge potential to 
the glue potential in full QCD by applying Eqs 
and 



(25) 



( |13[ ) improves the description of the chiral and 
(de)confinement phase transition with effective models. 



TABLE V. Pseudocritical temperatures for the crossover 
transition at ^/ = for different critical temperatures of the 
glue potential. They are determined by the peaks of the chiral 
susceptibility dAi^s/dT. 



[MeV] 


180 


210 


240 


270 


Tc [MeV] 


152 


158 


164 


171 



B. Dependence on the glue critical temperature 

In the last section we have showed that replacing the 
pure Yang-Mills Polyakov-loop potential by the quark- 
improved potential from functional calculations smoothes 
the phase transition significantly and brings it to closer 
agreement with lattice and functional calculations. Now, 
we will investigate the dependence of our results on the 
critical temperature of the glue potential T^^^'^. It is not 
yet uniquely defined so that we consider it as a parame- 
ter. According to the discussion in Sec. [TT| around Eq. ( 12 ) 
we expect it to be in the interval (180 - 270) MeV. For 
this investigation we use again the logarithmic parametri- 
sation of the Polyakov-loop potential and a medium 
sigma meson mass of TOq- = 500 MeV. For purpose of 
illustration we omit in the following figures to show the 
results with the Yang-Mills Polyakov-loop potential but 
choose the same range of the ordinates as in the previ- 
ous section to give an estimate of the reduction due to 
the application of the quark-improved Polyakov-loop po- 
tential. We show the results with both potentials in the 
same plots elsewhere [55] . 

Let us start again with the discu ssio n of the normalised 
pressure p/T^ 



perature of the Polyakov-loop potential is lowered. To- 
gether with the respective pseudocritical temperatures 
that are summarised in Table fV] the best agreement with 
lattice results shows the curve with a critical temperature 
of the glue potential of Ts'"° = 210 MeV. As discussed 
in the previous section the pseudocritical temperature of 
Tc = 158 McV for Tg'"' = 210 MeV is in accordance with 
the lattice resuh. But also the value for Tg""" = 180 MeV 
of Tc = 152 MeV is on the lower end of the transition re- 
gion of lattice calculations [HI] . As addressed already in 
the previous section including fluctuations to the model 
smoothes the phase transition additionally [3 and 
slightly increases the transition temperature as shows 
e.g. a comparison of Refs. and [S^. So we can ex- 
pect that taking into account fluctuations of quarks and 
mesons would decrease the transition temperature of the 
glue potential for which we find best agreement with the 
lattice calculation towards the 2-t-l flavour estimate of 
Ref. [6]. We will address the quantitative analysis in a 
future work. 

For the detailed investigation of the steepness of the 
transition in dependence of the critical temperature of 



We see in Fig. 14 that the pressure 
rises already at smaller relative temperatures and that 
it shows a slightly steeper increase when the critical tem- 



the improved glue potential we show in Fig. 15 the nor- 
malised trace anomaly or interaction measure. We see 
that using our adjustment of the Polyakov-loop potential 
towards the glue potential of full QCD the overall ampli- 
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FIG. 16. The Polyakov-loop expectation value $ as function 
of temperature at /i/ =0 with the adjustment to the glue 
potential for four different glue critical temperatures. The 
grey band is the continuum extrapolation of the lattice results 
of Ref. ^1] and the data points lattice results of Ref. [93j . 
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FIG. 17. The subtracted chiral condensate Ai,s as a function 
of temperature at /i/ = with the adjustment to the glue 
potential for four different glue critical temperatures. The 
grey band is the continuum extrapolation of the lattice results 
of Ref. [51] and the data points lattice results of Ref. [53] ■ 



tude is in good agreement with lattice calculations. Like 
anticipated in the discussion of the pressure the steep- 
ness of the transition increases when the critical tem- 
perature of the Polyakov-loop potential is lowered. Fur- 
thermore, the transition region broadens towards lower 
relative temperatures with decreasing transition temper- 
ature of the gauge potential. 



An important aspect that we see in Fig. 15 is that 
the slope of the trace anomaly becomes smaller and un- 
natural above for large critical temperatures of the 
Polyakov-loop potential. This behaviour has two rea- 
sons. First, the pseudocritical temperature is shifted to- 
wards the lower end of the transition region for larger 
transition temperatures of the glue potential, so that the 
transition region extends well above the pseudocritical 
temperature. So the emergence of the quark degrees 
of freedom holds on above Tc and avoids a decrease of 
the interaction measure. Second, the saturation scale of 
our adjustment between pure YM Polyakov-loop poten- 
tial and quark-improved glue potential, i.e., the upper 
limit of validity of Eq. (13 1, is closer to the transition 



scale for larger critical temperatures of the Polyakov-loop 
potential, see Table [Vj 

In Fig. [16] we show the dependence of the evolution of 
the Polyakov-loop expectation values on the critical tem- 
perature of the glue Polyakov-loop potential. The magni- 
tude of the transition increases when the transition tem- 
perature of the Polyakov-loop potential is lowered and 
the transition regions is shifted towards smaller relative 
temperatures. So the onset of deconfinement shifts to- 
wards the lower end of the transition region with decreas- 
ing critical temperature of the Polyakov-loop potential. 
This is a natural consequence, as decreasing the scale of 
the Polyakov-loop potential leads to a relative shift of the 
glue part towards a lower scale in the full theory. Over- 
all, the value of the Polyakov-loop variables shows a large 



sensitivity on the critical temperature of the Polyakov- 
loop potential but the general discrepancy to the lattice 
results is too large for any constraints. 

Due to the nontrivial coupling between Polyakov loop, 
quarks and mesons, the onset of deconfinement at lower 
temperatures with decreasing critical temperatures of the 
Polyakov-loop potential shows its impact also in the or- 
der parameter for chiral symmetry, the subtracted chiral 
condensate which is plotted in Fig. [T7| It deviates at rel- 
ative smaller temperatures from its vacuum value when 
the transition temperature of the Polyakov-loop potential 
is lowered. This effect gets enhanced due to the simulta- 
neous lowering of the pseudocritical temperature. Apart 
from that, the chiral condensate shows only a very mild 
dependence on the transition temperature of the glue po- 
tential. 



C. Dependence on the parametrisation 

The second uncertainty in the gauge sector of PNJL & 
PQM models besides the transition temperature of the 
Polyakov-loop potential is its form of parametrisation. 
We demonstrated in Figs.[8]and[9]that the different forms 
and parameter sets show a significant difference already 
in pure YM theory. Ideally, the Polyakov-loop potential 
of full QCD from first principal calculations should be 
applied. Since this is not yet directly applicable one can 
at least consider the uncertainty that comes along with 
the different parametrisations and parameter sets. 

For the comparison of the different parametrisations 
of the Polyakov-loop potential we choose a sigma me- 
son mass of 500 MeV and a transition temperature of the 
Polyakov-loop potential of 210 MeV. For better read- 
ability of the plots we omit in the following figures the 
results with the pure YM Polyakov-loop potential and 
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FIG. 18. Scaled pressure vjT^ as a function of temperature at 
fif — with the adjustment to the glue potential for three dif- 
ferent parametrisations and parameter sets of the Polyakov- 
loop potential. The orange band is the fit to the continuum 
extrapolation of the lattice results of Ref. [92) . 



choose the same range of the ordinates as in Sec.lTvllto 
give an estimate of the reduction due to the application of 
the quark-improved Polyakov-loop potential. The figures 
with the results with both potentials will be presented in 
Ref. IMl- 



The normalised pressure (Fig. 18) shows the steepest 



rise in the transition region for the polynomial parametri- 
sation with the parameters of Ref. [30] (Poly-I poten- 
tial). This is in accord with the description of pure 
Yang-Mills theory discussed in Sec. |IIIB| Overall, the 
difference in the results for the scaled pressure for the 
different parametrisations is not large around the transi- 
tion region. 

As seen before, the differences are more pronounced for 
the trace anomaly or interaction measure (Fig. 19). The 
slopes of the curves are similar for the Poly-I and loga- 
rithmic parametrisation but with a larger amplitude in 
the case of the Poly-I potential. This originates from 
the description of pure Yang-Mills theory as seen in 
Fig. [9] The normalised trace anomaly with the Poly- 
II parametrisation of Ref. [33] is nearly constant above 
the transition temperature. This is a feature of the offset 
of the maximum of the trace anomaly in Yang-Mills the- 
ory seen in Fig. |9] This offset is also seen in calculations 
of full QCD with the Yang-Mills Polyakov-loop potential 
[15] . The overall reduction of the amplitude in the tran- 
sition region for the improved glue potential leads to the 
behaviour seen in Fig. [19] 

The evolution of the Polyakov-loop variables shown in 
Fig. [20] reveals that the onset of the transition in the 
confined phase is the steepest with the logarithmic poten- 
tial. The Poly-II parametrisation leads to the smoothest 
evolution of the Polyakov loop below the pseudocritical 
temperature. Note that we normalised the Polyakov-loop 
expectation value shown in Fig. [20] for the calculations 
with the polynomial parametrisation such that it tends 



FIG. 19. Scaled trace anomaly (e — 3p) /T* as a function 
of temperature at ^/ = with the adjustment to the glue 
potential for three different parametrisations and parameter 
sets of the Polyakov-loop potential. The black band is the 
fit to the continuum extrapolation of the lattice results of 
Ref. [92] and the data points lattice results of the HotQCD 
collaboration '94'. 



towards unity at infinite temperature. Overall, the differ- 
ences in the Polyakov-loop variables due to the different 
parametrisations are so large that it should be considered 
before making quantitative statements. 

Nevertheless, the differences in the gauge sector due to 
the different parametrisations are considerably smaller 
in the evolution of the subtracted chiral condensate 
(Fig. 



21). The slightly stronger transition in the gauge 



sector for the Poly-I potential shifts the pseudocritical 
temperature to a smaller value compared to the other 
cases as seen in Table IVIl 
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FIG. 20. The Polyakov-loop expectation value $ as function 
of temperature at /// = with the adjustment to the glue 
potential for three different parametrisations and parameter 
sets of the Polyakov-loop potential. The grey band is the 
continuum extrapolation of the lattice results of Ref. [91] and 
the data points lattice results of Ref. [93] . 
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FIG. 21. The subtracted chiral condensate Ai,s as a function 
of temperature at /i/ = with the adjustment to the glue 
potential for three different parametrisations and parameter 
sets of the Polyakov-loop potential. The grey band is the 
continuum extrapolation of the lattice results of Ref. [ST] and 
the data points lattice results of Ref. [93], see also Ref. |95) . 
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FIG. 22. Scaled pressure p/T* as a function of temperature 
at /i/ = with the adjustment to the glue potential for three 
different masses of the sigma meson. The orange band is the 
fit to the continuum extrapolation of the lattice results of 
Ref. [H]. 



TABLE VI. Pseudocritical temperatures for the crossover 
transition at /i/ = for different parametrisations and param- 
eter sets of the Polyakov-loop potential. They are determined 
by the peaks of the chiral susceptibility dAi,s/dT. 





Poly-I 


Poly-II 


Log 


Te [MeV] 


144 


158 


158 



D. Dependence on the sigma meson mass 



Usually, the chiral partner of the pion is associated 
with the resonance /o(500) [ST]. Note however that in 
Ref. ^97) it was demonstrated that within an extended 
quark-meson model that includes vector and axial-vector 
mesons the resonance /o(1370) was identified as the non- 
strange scalar quarkonium state. Within our simple 
quark-meson model we consider to^ = (400 — 600) MeV 
to be a reasonable parameter range. 

To investigate the role of the mass of the scalar sigma 
meson the logarithmic parametrisation of the Polyakov- 
loop potential with a transition temperature of 210 MeV 
is used in the following. 

The normalised pressure as displayed in Fig. |22| shows 
that the steepness of the transition is relatively inde- 
pendent on the mass of the a-meson. The location 
of the pseudocritical temperature relative to the tran- 
sition region leads to differences in the evolution. For 
rria = 400 MeV the pseudocritical temperature is small 
(see Table VII) and at the lower end of the transition 
region so that the pressure just starts to rise at T^. The 
evolution for = 600 MeV is the other extreme. Here, 



TABLE VII. Pseudocritical temperatures for the crossover 
transition at fif =0 for different masses of the sigma me- 
son. They are determined by the peaks in the temperature 
derivatives of the subtracted condensate Ai,s . 



[MeV] 


400 


500 


600 


Tc [MoV] 


144 


158 


173 



the scaled pressure rises already significantly below 
since the pseudocritical temperature is just reached at 
the upper end of the transition region and at a relatively 
large absolute value as given in Table [VlTl 

The change of the pressure for an increasing mass of 
the sigma meson (Fig. 22 1 is opposite to the case of an 



increasing glue critical temperature (Fig. [14[). Therefore, 
to get agreement with the pressure of the lattice calcu- 
lation for a sigma meson mass of 400 MeV one would 
have to choose a critical temperature of the glue poten- 
tial of 180 MeV and for mo- — 600 MeV one would require 
"° — 250 MeV. So the lower the mass of the sigma me- 
son, the smaller the transition temperature of the glue 
potential has to be in order to reproduce the lattice re- 
sults. This is due to the fact that a decrease in the mass 
of the sigma meson lowers the scale of the chiral transi- 
tion as can be seen in Table VII and Ref. (321 ■ So the 



scale of the (de)confinement transition has to decrease as 
well. Note that the combination of m„ = 500 MeV and 
^ct"° = 210 MeV leads to a pseudocritical temperature 
that is equal to the one of lattice calculations. 

Figure [23] shows the results of the trace anomaly for dif- 
ferent sigma meson masses. The transition region broad- 
ens for a larger mass of the cr-meson. 

The subtracted condensate is shown in Fig. [24] The 
larger the mass of the sigma meson, the earlier the chiral 
condensate deviates from its vacuum expectation value. 
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FIG. 23. Scaled trace anomaly (e — 3p) /T* as a function 
of temperature at /i/ = with the adjustment to the glue 
potential for three different masses of the sigma meson. The 
black band is the fit to the continuum extrapolation of the 
lattice results of Ref. [92] and the data points lattice results 
of the HotQCD collaboration [94] . 
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FIG. 25. The Polyakov-loop expectation value $ as function 
of temperature at /i/ = with the adjustment to the glue 
potential for three different masses of the sigma meson. The 
grey band is the continuum extrapolation of the lattice results 
of Ref. [M] and the data points lattice results of Ref. [93] . 
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FIG. 24. The subtracted chiral condensate Ai_s as a function 
of temperature at /i/ = with the adjustment to the glue 
potential for three different masses of the sigma meson. The 
grey band is the continuum extrapolation of the lattice results 
of Ref. |91) and the data points lattice results of Ref. |93j . 



The temperature dependence of the chiral condensate be- 
comes smoother with larger tr-meson masses. 



The evolution of the Polyakov-loop variables (Fig. 25 ) 
shows that the transition regions shrinks when the mass 
of the sigma meson is lowered. Decreasing the glue crit- 
ical temperature broadens the transition region towards 
lower temperatures. The pseudocritical temperatures for 
Too- = 400 MeV is already smaller than the lattice result 
for rsi"<= = 210 MeV. 



In this work we have presented results for the non- 
perturbative Polyakov-loop potential in QCD in terms 
of the expectation value (^o)- The potential y[(^o)] 
is an order-parameter potential for the dcconfincmcnt 
transition. It is obtained from the first-principle QCD 
flows in Refs. [H [S] solely from the QCD two-point func- 
tions of the ghost and the gluon fields. In particular this 
computation takes into account the full backreaction of 
the quarks and the nonlinear coupling of the quark-gluon 
sector. So far, this important physics has been omitted 
in Polyakov-loop potentials used in PNJL/PQM models, 
where a Polyakov-potential is used which is generated 
solely from lattice Yang-Mills data for the Polyakov loop 
and the thermodynamic observables. 

In Sec|llj we have compared the glue Polyakov loop po- 
tential with that in pure Yang-Mills theory [THl [TS] [12] . 
We have shown that the two potentials are simply re- 
lated by an appropriate rescaling of the temperature, see 
Eq. ( 25 1 . Such an approach allows to mimic the effect 



of the quark backreaction on the gauge sector describing 
the confining dynamics. Whereas the fit of the param- 
eters of the Polyakov-loop potential to the lattice data 
contains only information about the minimum of the po- 
tential in the deconfined phase, our fit also adjusts the 
potential in the confined phase. In total this provides a 
simple method for improving effective low-energy models 
towards full QCD. 

In the second part of this work, we have applied this 
improvement in the 2+1 flavour PQM model. Our ad- 
justment of the temperature scale of the Polyakov-loop 
potentials leads to a softening of the temperature de- 
pendence of physical observables at temperatures close 
to the chiral and deconfinement crossover at vanishing 
chemical potential. As a consequence, the behaviour of 
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the order parameters and thermodynamic observables is 
in much better agreement with latest results from lattice 
QCD studies and first-principles functional approaches. 
Specifically, we find good agreement with the scaled pres- 
sure and the trace anomaly from lattice QCD simula- 
tions, even for the amplitude of the trace anomaly. 

In order to test the robustness of our results, we have 
analysed the dependence on the parametrisation and 
the transition temperature of the Polyakov-loop poten- 
tial as well as the cr-meson mass. We have found a 
parametrisation-dependence of the results, which can not 
be ignored in a quantitative analysis. To be more specific, 
we have varied the glue critical temperature and the mass 
of the scalar cr-meson and observed that a larger glue 
critical temperature has a similar effect on our results as 
a decrease of the mass of the cr-meson. In comparison 
to the computation without QCD-improved Polyakov- 
loop potential we find a qualitatively better agreement 
with the lattice data, see in particular the trace anomaly 
in Fig. [TT] For a sigma meson mass ma ~ 500 MeV 
and a transition temperature of the glue potential of 
~ 210 MeV this agreement is most pronounced. 

We have also discussed the next steps for systemati- 
cally improving PNJL/PQM models towards low-energy 
QCD. These steps include the determination of the in- 
put parameters in the chiral sector, i.e., mg., from the 
QCD flows. Furthermore, the model potentials can be 
upgraded by successively taking into account more infor- 



mation from the full glue potential. One possibility is e.g. 
to work-in the correct height of the barrier and the dif- 
ference of the potential at the minima in the metastable 
region. This has a direct impact on the steepness of both 
transitions. Finally the QCD potential obtained from 
the flow equation [H [5] can be used directly, see also 
Refs. ini Uni El] • Moreover, matter fluctuations can be 
included within an extension of the present work and 
Ref. [H]- 

In conclusion, the present work offers a simple and sys- 
tematic approach for improving the gauge sector of low- 
energy effective models towards QCD. The comparison 
of our results with that from lattice QCD confirms that 
this is important as it significantly affects the behaviour 
of physical observables in the crossover region. 
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